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Abstract 

Heavy flavours production in e-p DIS is studied at intermediate 
values of the transferred four-momentum square, under the assump- 
tion of boson-gluon- fusion mechanism dominance (no intrinsic heavy 
flavours contributions) . 

In this framework different expressions for the splitting functions in 
the gluon density evolution equation, with respect to the standard 
(DGLAP) ones, are explicitly derived. 
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Lepton-proton deeply inelastic scattering can be studied per- 
turbatively in the single-photon-exchange approximation neglecting 
electroweak corrections (Z exchange) when Aq CD < < Mf 
(q 2 = -Q 2 is the square of the transferred four momentum). 
The following expression for the differential cross section, with re- 
spect to Q 2 and the Bjorken scaling variable be obtained by 
Lorentz-invariance and current conservation arguments: 

2T ° 2 '2,,r/ ' ; (x, Q 2 ) + 2 (1 - y) F 2 (x, Q 2 )} (1) 
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a is the e.m fine structure constant (a ~ -^^j and y is the fraction 
of energy lost by the lepton, during the collision, in the proton rest 
frame. 

The proton electromagnetic structure functions F 12 (x, Q 2 ) satisfy 
the following factorization formulas [Q] 
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where p is the renormalization scale and pi is the density of the 
partonic species % in the proton with respect to the fraction of lon- 
gitudinal momentum carried by the parton itself (£), while e stands 
for the infrared cutoff used to regularize collinear divergences. 
The coefficients C[ 1 ^ (|, ^r, a s {p 2 )^j take into account the short- 
distance effects, they are infrared safe because of a Bloch-Nordsieck 
type compensation of soft divergences || and can be computed 
in perturbative QCD. 

The long-distance dynamics is enterely factorized in the parton den- 
sities pi (£,e,/i 2 ). 

The scale independence condition for the structure function F 2 

dF 2 (x,Q 2 ) _ Q 
d log p 2 

leads to the renormalization group equations for the quarks densi- 
ties, id est 
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where a s (/i 2 ) is the strong running coupling constant. 

The explicit expressions of the splitting functions Pff (z) and Pf g (z) 

are found to be IB] H 



P ff (z) = C 2 (F) 
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where in general [F (z)] + is a distribution defined as 



dz [F(z)] + <j>(z) 
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Charm production in e + p DIS has been extensively studied by ZEUS 
and HI collaborations at HERA electron-proton collider @ 
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In particular the charm contribution to the electromagnetic struc- 
ture functions of the proton is estimated from the measurements 
of the cross sections for inclusive D*^ production, given the value 
for the hadronization fraction of a quark charm in D* + by OPAL 
Collab. (e + e~ annihilation HIT 



D* 



0.222 ±0.014 ±0.014 



Experimental data can be explained in perturbative QCD by as- 
suming that the proton wave function does not possess any heavy 
flavour content; under this hypothesis heavy quarks production in 
e-p DIS is thought to be governed by boson-gluon-fusion (BGF) pair 
production according to the partonic subprocess (figure 1) 

7* 9 -> Qh Qh X 

That kind of experiments provides important informations on the 
gluonic structure of the proton; actually the most precise measure- 
ments of the gluon density in the proton available nowadays are 



coming from HERA experiments [Kj (T2[ . 
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The differential cross section for inclusive production of a pair of 
heavy quarks in virtual photon-proton scattering 

7* (q) p(P)^q h q h (M 2 ) X 

in pair invariant mass kinematics, enjoys the following factorization 
theorem 

daW (x, M 2 , Q 2 , m\) x 2 _ r 1 d£ 
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where represents the heavy flavour mass, M is the invariant 
mass of the pair while the sum is understood to be over all massless 
partonic species (gluon and light (anti) quarks). 

The coefficient functions a s (fJ 2 )^ include, as usual, 

the short-distance dynamics; they turn out to be infrared safe (com- 
pensation of soft divergences) and can be computed in perturbative 
QCD. 

The collinear divergences are regularized and factorized in the par- 
ton densities (£, e, /i 2 ). 
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Figure 1: BGF mechanism 

All physical quantities turn out to be scale independent, then 
from the condition 



d da( hh ) (x,M 2 ,Q 2 ,m 2 h 



d log// 2 dM 2 
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one may obtain the renormalization group equation for the density 
of gluons; at leading order it gives 
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The Born term of the lepton-proton cross sections for heavy flavours 
production in deeply inelastic scattering coincides with the lowest 
order contribution of the transition 

e g ^ e q h q h 

It comes from the partonic diagrams in figure 2: 




Figure 2: Born terms for e g — > e qh qh 



The calculations are straightforward and have already been per- 
formed in references [113] 

This contribution to the double-differential cross section can 
be written in the same way as in eq.(l): 



d 2 a 



'hh) 



2'KO 1 



dxdQ 2 xQ 4 



2xy 2 Ff h (x, Q 2 , m h ) + 2 (1 - y) F 2 hh (x, Q 2 , 



(8) 



5 



where (x, Q 2 , rrih) are the lowest order contributions of the heavy 
flavour h to the e.m structure functions of the proton. 
They can be expressed as follows 




while Qh is the heavy quark electric charge (in e units). 
The aim of this paper is to derive the explicit form of the splitting 
functions P g f (z) and P gg (z) which appear in the evolution equation 
for the gluon density, directly from the study of heavy flavours pro- 
duction in e-p DIS. 

Then these results are compared with the corresponding quantities 
obtained by Altarelli and Parisi in a different framework || [[J. 
For this purpose we start by considering the next-to-leading order 
(NLO), quark initiated corrections to the Born term, namely the 
lowest order contributions to 

e qj -> e q f q h q h 
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where qj is here any light quark. They are given by the diagrams 
shown in figures 3 and 4 (all the calculations in the present paper 
are made in Feynman Gauge). 

The corresponding contribution to the differential electron-proton 
' /2<T 2 cross section is expressed, in the framework of free parton 



dxdQ 

model, as a convolution over the free (light)quark density 



dxdQ 2 J x £ dzdQ 2 

sents the "partonic" Bjorken 
is the corresponding lepton-parton cross section. 



where z — § represents the "partonic" Bjorken variable and — - — 





Figure 3: NLO quark-initiated corrections; the e.m 
interaction involves the light quark 



Figure 4: NLO quark-initiated corrections with the 
virtual photon probing a heavy flavour 
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One can easily realize that the collinear divergent part of dxd q 2 
comes from the last two diagrams in figure 4; it can be written in 
the same form as the Born contribution 



dxdQ 2 Jo £ P dzdQ 2 



P 9 (0 ( 9 ) 



simply by subsituting the "free" gluon density p g (£) with its next- 
to-leading order correction |L5f 



where is the infrared cutoff used to regularize the light quark- 
gluon collinear singularities (we choose Q as factorization/renormali- 
zation scale). 

Therefore we can regularize these NLO collinear divergences and re- 
absorb them in the Born term through the following redefinition of 
the gluon density in eq.(9) 

mo ^ p 9 (0 + (e, |j) (io) 

The splitting function P g f (z) is given by 

P gf (z) = C 2 (F)z = ^z (11) 

It should be noticed that this expression vanishes linearly in the limit 
z — > and that it is different from the following formula obtained 
in references IHI [El 



P gf (z) = C 2 (F) 



z 

To compute now the explicit expression for the gluon-gluon splitting 
function P gg (z), the next-to-leading order, gluon initiated correc- 
tions, due to real gluon emission, have been derived in the collinear 
limit. 

Essentially the following partonic process has to be studied (at low- 
est order): 

e g -> e g q h q h 

There are in total 8 diagrams to consider, they are shown in figures 
5 and 6: 
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Figure 5: NLO gluon-initiated corrections involving 
only quark-gluon vertices 
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Figure 6: NLO gluon-initiated corrections containing 
the three gluon vertex 



The contribution of those diagrams to the electron-proton dxd q 2 
cross section can be written, in the free parton model, as a convo- 
lution over the free gluon density: 



Ihh] 



dxdQ 2 
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where dzd q 2 is the corresponding lepton-gluon cross section. 
The leading term in the gluon-gluon collinear limit comes from the 
last 2 diagrams in Figure 6; furthermore it can be written in the 
same form as the Born contribution if one replaces the free gluon 
density p g (£) in eq.(9) with the next-to- leading order correction 

Q 2 \ a-, (Q 2 \ f ldX „ (£ 
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where the explicit expression for the splitting function P gg (z), valid 
for z < 1, is given by |T5] 

z 



P gg (z) = 2C 2 (A) 



1 - z 



+ z(l-z) 



, C 2 (A) = N 



(12) 



Notice that eq.(12) differs from the standard expression derived in 
references || @, id est 

2C 2 (A) \-^- + z {l-z) + — 1 (13) 



p 99 w 
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Moreover we point out that P gg {z) in eq.(12) vanishes linearly as 
z — > and that it is not invariant under z 1 — z exchange contrary 
to the previous expression. 

In addition both P g f (z) in eq.(ll) and P gg (z) in eq.(12) can be ob- 
tained from the corresponding Altarelli-Parisi results by neglecting 
all terms that do not vanish as z — > 0. 

Therefore the NLO gluon-gluon collinear divergences can be regu- 
larized and then reabsorbed in the following redefinition of the gluon 
density in the "Born level" term (eq.(9)) 

p 9 (0 =► p 9 (o + s P f (e, |I) (14) 

Finally 5pf> (f, turns out to be divergent in view of a soft sin- 
gularity in the limit £ — > 1 (gluon emitted with arbitrarily small 
energy) . 

It is a peculiar feature of gauge theories due to the fact that gauge 
bosons are massless. 

If one examines the NLO interference terms of the virtual corrections 
with the Born contributions for e g —>■ e qh qh transition, related to 
collinear divergences, soft singularities are found, which compensate 
the analogous divergences in the real emission contributions. 
This mechanism is similar to the Bloch-Nordsieck compensation well 
known in QED gj. 

The mentioned virtual corrections are given by the six diagrams in 
figures 7 and 8 (three more diagrams are obtained from those in 
figure 7 by reversing the heavy quark line). 

The loops in figure 7 involve only massless states; they are indeed 
scaleless integrals which vanish if appropriately defined in dimen- 
sional regularization 0. 

Only three diagrams are left (figure 8); their interference terms with 



the Born contributions have been included in the calculation [13 . 
In conclusion the overall next-to- leading order correction to d ^ d g 2 , 
for the heavy flavour h production (via boson-gluon-fusion) , com- 
puted in the collinear limit, can be reabsorbed in the Born contribu- 
tion by means of the following redefinition of the free gluon density 
in eq.(9): 

p a (0 =► p 9 (0 + 5 Pg { e, %) (is) 



Q 
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with 



Q' 1 



q 2 \ r'dt 



The splitting functions P g f (z) and P gg (z) are given by 

1 



(16) 



P gf (z) = C 2 (F) z , P gg (z) = 2C 2 (A) 



'l-z) 



The sum is understood to be over all light (anti) quarks. 




Figure 7: Gluon self energy contributions of order a s 



Therefore the effective gluon density involved in heavy flavours 
production in e-p DIS, at intermediate values of the factorization 
scale, follows a different evolution equation with respect to usual 
DGLAP eqs.; needless to say this fact affects the (anti)quarks den- 
sities too as soon as the evolution equations are coupled. 



12 



This feature is due to the assumption on the dominance of boson- 
gluon-fusion mechanism; it is indeed a consequence of the explicit 
flavour symmetry breaking (intrinsic heavy flavours contributions 
are suppressed) and it is expected to be verified at moderate values 
of the scale /i with respect to m h . 

For n values much higher mass effects are negligible and the heavy 
flavour h behaves as a light (massless) flavour; the intrinsic contri- 
butions to the cross sections can no longer be neglected. In that 
case the gluon density follows again the standard DGLAP evolution 
equation so as the density for h quarks. 




Figure 8: e g — > e q^ qn virtual corrections 
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